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We introduce a systematic method for extracting multivariable universal scaling functions and 
critical exponents from data. We exemplify our insights by analyzing simulations of avalanches in an 
interface using simulations from a driven quenched Kardar-Parisi-Zhang (qKPZ) equation. We fully 
characterize the spatial structure of these avalanches- we report universal scaling functions for size, 
height and width distributions, and also local front heights. Furthermore, we resolve a problem that 
arises in many imaging experiments of crackling noise and avalanche dynamics, where the observed 
distributions are strongly distorted by a limited field of view. Through artificially windowed data, 
we show these distributions and their multivariable scaling functions may be written in terms of 
two control parameters, the window size and the characteristic length scale of the dynamics. For 
the entire system and the windowed distributions we develop accurate parameterizations for the 
universal scaling functions, including corrections to scaling and systematic error bars, facilitated by 
a novel software environment Sloppy Scaling. 



PACS numbers: 64.60.De, 89.75.Da, 75.60.Ej, 05.10.Cc 
I. INTRODUCTION 

Systems that have crackling noise and avalanches ex- 
hibit scale invariance and power laws, which point to 
the notion of underlying universality [1 j. These sys- 
tems include many of the best-studied examples of non- 
equilibrium critical phenomena, and much progress has 
been made in a renormalization group context [2j [3] . The 
renormalization group implies that the long length and 
time behavior near critical points is governed by universal 
exponents and scaling functions. However, the predictive 
power of these theoretical studies has hiterto been under- 
utilized; the primary focus of experiments and numerical 
simulations has been on precise estimates of critical ex- 
ponents, rather than on the universal joint predictions of 
properties involving several control parameters and/or 
measured quantities. 

A wide variety of materials and natural systems have 
been studied in the context of non-equilibrium critical 
phenomena. Many of these systems exhibit avalanches 
which have power law size distributions. These include 
Barkhausen noise in ferromagnets [4-8 , fluid imbibition 
into porous media, flux-line depinning [3] I9HT3], and 
martensitic transformations [HI [15], to name a few. In 
the first three of these systems, avalanches are the result 
of the jerky motion of an interface (domain wall, fluid 
front, flux-line) in a disordered environment, and can be 
described by the same family of front-propagation mod- 
els. 



In this manuscript, we study the spatial structures of 
avalanches in a front-propagation model in two dimen- 
sions, developing tools and methods needed for system- 
atic study and extraction of these multiparameter univer- 
sal scaling functions. To illustrate the utility and impor- 
tance of these functions, we apply them systematically to 
a practical experimental problem - the size distributions 
of avalanches seen through a viewing window. This prob- 
lem illustrates (a) the complexity and sophistication of 
the different emergent size distributions, (b) the relation- 
ships amongst the probability distributions of heights, 
widths, and sizes and their utility in generating predic- 
tions for windowed avalanches, and (c) the use of func- 
tional forms and least-squares fits to analyze and report 
on these multiparameter scaling functions. 

Imaging experiments have been used in recent years to 
study a wide variety of systems exhibiting crackling noise 
or similar dynamics. Barkhausen noise is measured mak- 
ing use of the magneto-optical Kerr effect [5 , 6 , 16 , allow- 
ing one to examine the domain wall motion in 2D thin 
films in real-time. In experiments on superconducting 
vortices [17], a magneto-optical (MO) setup is also used. 
In experiments of fracture [18], fluid imbibition [19] and 
granular systems [20] [21] , the dynamics are also followed 
with a camera. 

These visualization experiments provide an unusual 
opportunity: we now can study the universal proper- 
ties of the spatial morphology - various distributions of 
heights, widths, angles, local heights, etc of either the 
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avalanches or the fronts. However, the measurements 
of these properties are often distorted by a limited field 
of view. We hereby take this problem and develop the 
scaling theory for the universal functions needed to char- 
acterize the results of a generic imaging experiment - the 
distribution of avalanche sizes seen through windows. 

The limited field of view in experiments distort the 
size distributions of avalanches, and cause difficulties in 
characterizing the critical exponents. Naturally, there is 
a bias towards small avalanches; large ones are cut off 
by the boundaries of the window. It can also distort 
the size distribution if pieces of large avalanches cut off 
by window boundaries are counted as small ones. Ex- 
periments have taken a variety of approaches to deal 
with such windowing effects: "Laser reflect ometry" [16] 
on Barkhausen noise uses the magneto-optical Kerr ef- 
fect, but lacking spatial resolution, lumps fragments and 
avalanches together; meanwhile, other optical Kerr ex- 
periments have shown [6 that the effective size expo- 
nent r for this lumped distribution depends strongly 
on the window width. Work by Kim et al. [22] report 
quite striking distributions but do not specify whether 
their data includes avalanches that touch the bound- 
ary. In experiments on superconducting vortices, or mag- 
netic flux avalanches [17] avalanches exceeding a certain 
size are discarded. In fluid imbibition [19] , the edges 
of the system are purposely left out to avoid any dis- 
tortion produced by side walls. In granular systems, 
where avalanche dynamics in piles of rice are studied 
with real-time reconstruction [20] [21], and fracture ex- 
periments [T8] , where the dynamics are followed with a 
high-speed camera, boundary effects are not considered 
but may also be important. 

We will show comprehensively how to analyze all of the 
size data lying within a window, and how to use the differ- 
ent classes of avalanches to get independent measures of 
various critical exponents. Indeed, window-width finite- 
size effects need not be avoided, but properly treated may 
provide additional measures of the critical exponents and 
the spatial structure. 

Characterizing spatial structures of avalanches must 
go far beyond the traditional focus on critical exponents. 
Many experiments report power laws, however through 
this study we emphasize that one can make predictions 
about both power laws and scaling shapes from data, as 
has been demonstrated in a previous study of avalanche 
temporal shapes [23] . Indeed, traditional scaling collapse 
methods fail for functions of more than two variables. 
To optimally extract the behavior and estimate errors, 
we need to do simultaneous analysis of many different 
properties. We thus introduce a software environment, 
Sloppy Scaling [24 , which facilitates the exploration and 
development of simultaneous fits of multiple data sets 
with parameterized forms of universal scaling functions. 
With this approach we are taking the first steps towards 
the use of scaling methods as a practical engineering tool. 



II. SUMMARY OF KEY RESULTS 

Since our theme is multifaceted, readers may be inter- 
ested in focusing on different aspects of this work. In 
this section we present an overview of key results in this 
paper and a summary of their corresponding sections. 

1. Universal spatial structures of avalanches 
in directed percolation depinning. We provide a 
substantive analysis of the universal spatial morphol- 
ogy of avalanches in the quenched KPZ (qKPZ) model 
in 1+1 dimensions (the model is discussed in section 
III). Figure [TJ a) shows avalanches in a typical simula- 
tion of this model. Analogous to magnetic systems, we 
have added a "demagnetization factor" k that param- 
eterizes the restoring force, which allows us to access 
many metastable configurations of the front near the de- 
pinning transition, and controls the typical width of an 
avalanche, ~ k~ Vh . Figure [2] shows examples of these 
resulting avalanches from simulations of various k. In sec- 



tion IV we thoroughly examine the spatial structure of 
the avalanches, including sizes s (the total area covered 
between pinned fronts, which would correspond to the 
total magnetization change in magnets, or the avalanche 
size of a rice pile) , and also widths w and heights h (which 
measure the length of an avalanche in directions perpen- 
dicular and parallel to the direction of the motion of the 
front; this corresponds to studying the shapes of the mag- 
netic domains or flux lines). We examine and fit the dis- 



tributions of these sizes, heights and widths in section IV 
2. Avalanches in Windows. As mentioned in the 
Introduction, in many imaging experiments the limited 
field of view distorts the avalanche size distribution. This 
is illustrated in Fig.[TJb) and (c). In Fig.[TJb) most of the 
avalanches are cut off by the left and right boundaries- if 
one were to count the area of these avalanches, we would 
count large ones as smaller ones. The resulting size dis- 
tribution inside the window (dashed line of Figure [TJc)) 
has a very different power law and shape compared to 
the full system distribution (solid line of Figure [ljc)) . 

In section |V| we show in detail how avalanches which 
cross boundaries exhibit distinctly different size distribu- 
tions and critical exponents. For simplicity, we consider 
a strip geometry where only the left /right boundaries 
are relevant (since avalanches in our system are flat and 
wide, few touch the top or bottom). We can therefore 
separate avalanches into four different categories, 00, 01, 
10, and 11, depending on whether avalanches touch (1) 
or do not touch (0) the left or right boundaries: internal 
avalanches (00) (Figure [TJd)), split avalanches (10 and 
01) (Figure file) and (f)), and spanning avalanches (11) 
(Figure [ljg)). The internal (00) avalanches share the 
same power law as that of the full distribution, with a 
cutoff controlled by both the system size and the win- 
dow size W. The split (10 or 01) avalanches have a mod- 
ified power law- a smaller exponent as larger avalanches 
are counted as smaller halves. (See Figure [T])(c)). The 
spanning avalanches also exhibit a smaller exponent, al- 
though this is not obvious in the shape of its scaling 
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FIG. 1: Windowed distributions (color online) (a) The full system of avalanches. This figure shows a qKPZ simulation 
with avalanches, (b) A limited field of view: we can only see part of the system. The avalanches inside the artificial viewing 
window are brightly colored, and those outside are washed out. Notice that the avalanches within the window are cut off at top 
and bottom, and (more importantly for short, wide avalanches with roughness exponent £ < 1) on the two sides, (c) The size 
distributions for the different types of avalanches: (d) internal 00 avalanches, (e) split 10 avalanches, (f) split 01 avalanches, 
and (g) spanning 11 avalanches 



function. We can see that it has both an outer cutoff 
due to Lk and an inner cutoff due to the window size 
W, since avalanches must be large enough to span the 
window (purple dash-dotted line in Figure [TJc)). The 
internal and spanning avalanches also have distinct uni- 
versal scaling functions with a cutoff controlled predomi- 
nately by the window size for windows comparable to or 
smaller than the size of In Section |V| and Section 
VI] we give a thorough analysis of these modified power 
laws, the different scaling shapes, and the results of fits 
to data. 

Having established a sophisticated method of analyz- 
ing both experimental and simulation data, we can utilize 
this analysis to enhance the collection of data in visual- 
ization experiments. Section |VII| has some suggestions 
for how to collect data and simultaneously analyze the 
scaling behavior of different magnifications, and extract 
multiple exponents. 

3. Functional forms A main emphasis of our work 
is that we fit an entire functional form instead of power 



laws [25], this includes the shape of the scaling function, 
and analytic and singular corrections to scaling. The ben- 
efit of approaching a scaling problem this way is that it 
allows us to account for both universal and non- universal 
effects in a consistent way. Writing down the functional 
form that is given by the data for a certain universality 
class will also be useful for identifying and characteriz- 
ing other systems that are thought to belong to the same 
universality class. 

We have found that, to analyze the windowed distri- 
butions, we need to first thoroughly examine the spatial 
structure of avalanches for the full system. In partic- 
ular, in order to analyze the avalanche pieces left inside 
the window, we need to define height and width distribu- 
tions and also joint distributions of sizes and widths. Sec- 
tion |IV] discusses these and also the results of fits for such 
distributions for the qKPZ model we are studying. To fo- 
cus on the scaling region, and minimize lattice effects, we 
will discuss these distributions in terms of fractional area 
distributions, A(S) oc SP(S), the average fraction of the 
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system that a given size takes up. 

A remarkable result is that size, height, and width dis- 
tributions can be fitted with a nearly identical functional 
form. For example the size distribution is: 



A(S\L k ) = (S/Ll +( ) 2 -TAs k (S/Ll + <)/S 



(1) 



ZsSi s ))/S. 



Here we have the shorthand S k = S/L 1 ^^ . The width 
A(w\L k ) and height A{h\L k ) distributions are similar, 
with the form of ±Y a A(y), where X = {S,h,w} and 
Y = {S/ ' L 1 ^^ , h/L^w/Lk}. A(y) is identical in form to 
the one quoted above in Eq[T] 

Fitting the size, height and width distributions at once, 
we can extract multiple exponents- not only the com- 
monly measured size distribution exponent r, but also 
the exponent v k which, as mentioned in point 1, relates 
the typical width of an avalanche to fc, and also the rough- 
ness exponent £. The roughness exponent £, which mea- 
sures the fluctuations of the interface, is typically quoted 
for front propagation models as a means of character- 
izing the universality class. (In our analysis, we have 
found a range for the roughness exponent ( (from 0.62 
to 0.72) that differs from the literature value 0.63 [13]; 
we discuss this in Appendix |b|) One must note that the 
functional forms we choose are a practical tool to sum- 
marize existing information. While they may be inspired 
by analytical calculations, and chosen to be consistent 
with known asymptotics, they should be trusted only in 
the ranges over which they have been measured. 

4. Mult i variable scaling problems. The size dis- 
tribution in Eq. [I] has a scaling form with one scaling 
variable. However, in this paper we will consider many 
scaling forms with more than one variable, such as a joint 
size and width distribution (IV), or the windowed distri- 
butions (section |V| and VI). In these cases, two or more 
scaling variables are important for describing the shape 
of the distribution (as seen in Figure [TJc)). For example, 
the general form for the 11 windowed distributions is: 



AMWLk) 



1 (Sk\ 

s \Wk) 



(2-r)(l+C)/C 



An(s k ,W k ). (2) 



Here the scaling functions become distributions with two 
scaling variables, the rescaled size s k and rescaled window 
width W k . The traditional "scaling collapse" methods 
become problematic when multiple scaling variables are 
simultaneously important; this has hitherto retarded the 
effective study and use of these powerful universal joint 
distributions. 

We present a systematic method for analyzing scaling 
problems with multiple control variables. In our prob- 
lem, the two control parameters are the demagnetization 
factor /c, and the window width W. We will show that 
the interplay between k and W is important for deter- 
mining the shape of the avalanche size distributions. In 
particular, we can write scaling functions for the distribu- 
tions of the avalanche pieces in terms of these two scaling 



variables s k = s/L 1 ^^ and W k = W/L k (as seen in Fig- 
ure [TJc)). For example, the scaling function for the 11 
distribution is: 



Au(s k ,W k ) = exp(-(T„ + U u sl /2 + Znsl 11 



• Cn 



Sk 



-ml 



) (3) 



Section VI discusses in detail the functions Aoo(s k , W k ), 
Aio(s k , W k ), and A\\(s k ,W k ) that we consider. One 
may also find the results of fits (figures and tables of 
parameters) in section [V[ 

5. Sloppy Scaling The analysis in this paper is done 
in the software environment Sloppy Scaling and includes 
a Bayesian analysis of systematic error bars, which are 
explained in Appendix |A4| Sloppy Scaling allows us to 
fit data without collapses, which as mentioned above is 
problematic when there is more than one scaling vari- 
able involved. Included in the software setup are auto- 
matic fits of data to theory functions with nonlinear least- 
squares and ease of visualizing results. This software 
may be applied to many different multivariable scaling 
problems, making full use of universality and the pre- 
dictions of the renormalization group. All of the fits in 
this paper and their corresponding figures (including axis 
labels) were generated directly and automatically using 
Sloppy Scaling. 



III. MODEL 

We use a model for imbibition fronts to produce 
avalanches, it is extensively studied, with well-established 
critical exponents. 

We simulate crackling noise using a quen ched KPZ 
model in 1+1 dimensions (see Appendix |A 1| for details 
on the implementation) [10], with dynamics given by: 



dh(x, t) 
dt 



F - k(h) + 7V 2 /i + \(\7h) 2 + ri(x, h) (4) 



where h(x,t) is the height of the front, F the driving force 
increasing quasistatically, linear and non-linear terms for 
the KPZ model controlled by the parameters 7 and A 
respectively, and 77 Gaussian quenched disorder. In the 
spirit of magnetic avalanche models [26 , we have added a 
term analogous to a demagnetization field —k(h), which 
allows us to have pinned fronts at many metastable con- 
figurations. Models like these have been simulated [9 
mostly near the depinning transition (and with k = 0). 
With k 7^ 0, we define the area between each pinned 
front as an avalanche of size S. Avalanches produced by 
this model are thought to belong to the directed perco- 
lation depinning (DPD) universality class [I3j[27]. The 
avalanches are self-afflne, long and wide, with £ < 1. The 
roughness exponent ( characterizes the ruggedness of the 
front, and also governs the scaling of avalanche heights h 
with widths w, h (xw^ [9j|28]. In our model, k controls 
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law times a universal scaling function. For example: 



FIG. 2: Scale Invariance in Crackling Noise, (color 
online) Three simulations of crackling noise, with different 
"demagnetization fields", k. (a) k = 10 -5 (b) k = 10 -6 
(c) k — 10 -7 (k is what controls the typical avalanche sizes 
in a given system, giving a characteristic width Lk.) Larger 
demagnetization fields stop avalanches more strongly, hence 
large k corresponds to smaller avalanches. The colored re- 
gions represent avalanches. The fronts are moving from bot- 
tom to top. Notice that the two simulations are statistically 
similar to one another apart from a rescaling of heights and 
widths. Note that most of the area is covered by the largest 
avalanches. 



the typical size of the avalanche (Figure [2j, the larger the 
k the smaller the typical size of an avalanche. We define 
the characteristic width of an avalanche in the full, un- 
windowed system to be L^ = k~ Uk . We will write all the 
scaling forms in terms of the length rather than the 
demagnetizing factor k directly to emphasize the analo- 
gies to finite-size scaling, as we are also studying finite 
window sizes compared to the size of 



IV. AVALANCHE SHAPES AND 
DISTRIBUTIONS OF THE FULL SYSTEM 

In this section we will introduce various avalanche spa- 
tial distributions and their scaling forms. These forms 
will allow us to motivate the windowed distributions, and 
also serve as an example where traditional collapses may 
lead to questionable results. 



A. Area-weighted size distributions 

Traditionally, to describe an avalanche size distribu- 
tion, we write the probability distribution as a power 



P(S\L k ) = S- T P(S/Ll +c ). 



(5) 



Lk is the characteristic width of an avalanche, and h ~ 
the typical largest height, and therefore the appropriate 
scaling variable to describe the area of an avalanche is 
S k = S/L 1 ^. 

When studying the spatial properties as in our case, 
the probability distribution P(S) (Eq. [5| is not the best 
choice, as its normalization is highly affected by non- 
universal effects at the lattice spacing (Figure [3]). In fact, 
for r > 1, the normalization integral 
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P(S\L k )dS. 



S~ T dS 



,2(1- 



L 



(l-r)(l+C) 



(6) 



diverges at its lower (ultraviolet) limit a — >• 0, but not for 
L/e — )> oo. Although we could study scaling functions that 
include a lattice cutoff, it is more interesting to focus on 
the large (infrared) avalanche cutoff, which depends on 
L/e in a universal way. To this end, it is more appropriate 
to make use of the first moment of P(S), and to consider 
the area-weighted size distribution A(S) 



A(S)~S>P(S). 



(7) 



A(S)dS has a natural physical interpretation: it is the 
fraction of the full system area covered by avalanches 
with sizes between S and S + dS. Its scaling form is 
thus: 



A(S\L k ) = L k T - 2){1+ ^S^A S k{S/L\ +i ) 



= {S/L\ + ^Ask{S/L\^)/S 
= Sl~ T Ask(Sk)/ S. 



(8) 



We use the subscripts Sk to distinguish the scaling of 
the size distributions governed by k from those governed 
by other control variables. The power of Lk we pull out 
of the scaling function is needed to normalize A(S) to 
one, since A(S) is sensitive to the long-distance cutoff. 
In particular, 
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(a) Size distribution (probability) 



(b) Area- weighted size distribution 



FIG. 3: A(S) vs P(S). (color online) Here one can see the difference between the more traditional P(S) and the area- weighted 
A(S). (a) P(S): most of the area under the curve is from small avalanches, where non-universal lattice effects are important, 
(b) A(S): the normalization is dominated by large avalanches, avoiding the lattice effects, so that we can focus instead on the 
dependence on the large scale cutoffs- L k and, in later sections the window size. The data here is from qKPZ simulations of 
different k with different simulation size L. The lines in (b) are a result of a joint fit with the maximum height and width 
distributions. The fitting function and fitting parameters are shown in Table [I] 
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where the last integral converges for r < 2 [43] and be- 
comes small as L k becomes large. Notice that the nor- 
malization of a power law must either diverge at the lower 
or at the upper limit. By studying avalanches weighted 
by their first moment, the normalization depends explic- 
itly on L k , the infrared cutoff. This is the regime we 
are mainly interested in, since we would like to study 
the finite size effects imposed by both k and the win- 
dow size W. Alternatively, as done in [29 , one may also 

/C 2 \ 

define a scale S m = redefine the sizes as S / S m 

and a corresponding p(S/S m ) which has normalization 
J °° dssp(s) = 1. Here p(s) is universal, but is also 
not a probability distribution in the conventional sense. 
This definition has an effect equivalent to what we do 
here: to make the function universal and insensitive to 
non-universal lattice effects on normalization. Namely, 
our definitions are related in the following way: A(S) = 
^rp(S/S m ) and A Sk (S/S m ) = (S / S m ) 2 p(S / S m ). Here 

we use Sm = S k = S/L 1 ^, which is consistent with 
their definition of S m up to a constant factor. We pre- 



fer to focus on the more directly interpretable fractional 
area distribution A(S). 

Furthermore, notice that we have been unorthodox in 
writing the scaling form ([8| for A(S\L k ) with a power of 
both L k and of S outside the scaling function. Normally 
one factors out a single variable from the scaling function. 
For example, one could in principle write 

A{S\L k ) = ^B Sk (S k ). (9) 

In this form, A(S\L k ) dS is invariant under rescaling, and 
also clearly preserves normalization. However, B{S k ) = 
S 2 ~ T A(S k ) vanishes as S k — )► 0, so this form of the 
scaling function disguises the power law behavior of the 
avalanche size distribution. By choosing the form ^4^^ 
which is defined to be finite and non-zero as S —> 0, we 
make manifest both the avalanche size dependence and 
the system size dependence. 
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B. Maximum height and width distributions 

In addition to the size distributions, we can also study 
the avalanche height and width distributions. We define 
height along the direction of front propagation, measur- 
ing the maximum height of an avalanche, and width per- 
pendicular to heights, measuring the maximum width of 
an avalanche. 

How do the height and width distributions scale? An 
avalanche of height h has size S ~ hw ~ hh 1 ^ = 
ft( 1 +C)/C ? go the system area A(h\L k )dh covered by 
avalanches with heights between h and h + dh scales as 

A(h\L k )dh ~ A(S\L k )dS 

= 4--2)(iK) s i-r^ fc(s/i i+c )dS 

= L^^+O^d-rJd+O/C Ahk (h/Ll) dS 

(10) 

and since dS/dh ~ h 1 ^ , 

= (v4) (2_T)(1+0/< ^*(V4)A 

= /i (2-r)(l +C )/C Afe(/lfc)//l (n) 

where = ft/i^ is the appropriate scaling variable to 
describe the avalanche height. 

Similarly, we can write the scaling form for the width 
distributions as: 

A(w\L k ) = (w/L k Y 2 - T ^ 1+ <U wk (w/L k )/w 

= w {2-T ){ i+Q Awk{wk)/w (12) 

where w k = w/L k is the scaling variable for widths. 

Notice the pattern in all these formulas: A(X\L) is 
a function with A(Y) of scaling variables (combinations 
of X and L invariant under the renormalization group 
rescaling) multiplying a power of the scaling variable, di- 
vided by the independent variable X, making A{X\L)dX 
invariant under rescaling, allowing it to be universal. 



C. Joint fit of size, height and width distributions 



are listed in Table [I]). The first exponential term is the 
universal scaling function, while the second accounts for 
non-universal analytic corrections at small x due to lat- 
tice effects [30]. (See Appendix A 5). 



In Figs. [3(b)] [4(a)] [4(b)] we show the size, height and 
width distributions for our qKPZ simulations at various 
fc, and at different simulation sizes L. The curves are 
theoretical fits using a functional form of the type 

A(Y) =exp(U x Y 1/2 - Z x Y 6 *)exp(Af/x + A x 2 /x 2 ). 

(13) 

where x = {S,h,w}, and Y = {Sk,hk,Wk}, and 
U x , Z x , S x , Af , A% are (5x3) fitting parameters (results 



One may ask why we choose this particular scaling 
form. When fitting data to a function, there are many 
parameterizations we could use to describe the data. 
This form is motivated from a functional renormaliza- 
tion group expansion by Le Doussal and Wiese for static 
avalanche size distributions in a linear model [31 j. Our 
model differs in that there is a nonlinear KPZ term lead- 
ing to anisotropy, so that our avalanches belong to a dif- 
ferent universality class, the DPD universality class. 

Le Doussal and Wiese find the avalanche size distribu- 
tion for the linear model, for all static universality classes 
(random-bond, random-field, and random-periodic), to 
leading order in d = 4 — e (where d is the dimension of 
the interface), to be: 



P(S) oc S~ T exp(C(5/5 m ) 1 / 2 - ^(S/S m ) 5 ) (14) 



Here their scaling form includes the large scale cutoff 
S m . Le Doussal and Wiese claim that their results for 
both static and dynamic avalanches agree up to one loop 
for systems with £ < 1. Static avalanches are separated 
by equilibrium configurations (or ground states), and dy- 
namic avalanches are connected through a sequence of 
metastable states. Our avalanches result from a moving 
interface near the depinning transition, so they belong 
to a dynamic universality class. We thus use Eq. [14] as 

However, 



inspiration for the scaling form of Eq. 13l However, one 
may note that there is no theoretical basis that it should 
work, since it is from a distinctly different universality 
class. 

One cannot determine the values of v k and £ indepen- 
dently, if we fit the size, height, or width distributions 
with Eqns.[8j [TT] and [12] separately. For example in the 
size distribution we can only determine the combination 
v k (\ + C)- We determine the three critical exponents r, 
and C, by jointly fitting the size, height and width 
distributions. 

The results of our fits are reported in Tableland shown 
in Figures 3(b) and [4] In particular, we find ( = 0.62 ± 
0.02 which is close to the highly-precise value of £ = 0.63 
found in the literature [T2l l28l I32H34] . We note that the 
parameters 5 X , which in principle control the aymptotic 
decay of the scaling function, are estimated here from a 
fit to the entire distribution. The quoted errors do not 
represent a confidence on the asymptotic decay - merely 
a confidence in the predictions over the range where data 
has been fit. 
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FIG. 4: Area- weighted avalanche distributions for maximum heights and widths (color online) Area- weighted 
avalanche distributions of (a) maximum heights, and (b) maximum widths for qKPZ simulations at different k, and simu- 
lation sizes L (dots are binned data). The critical exponents were jointly fit with the size distributions A(S\L k ) of Fig. [3] 
using the scaling forms of Eqs. |12[ and |13| The best fit values for the critical exponents, parameters for universal scaling 
functions, and non- universal corrections are given in Table [I] 
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FIG. 5: Size distribution A(S\Lk) scaling collapse (color 
online) We collapse the size distribution data using the uni- 
versal exponents of r = 1.24, Vk — 0.45, and £ = 0.62, the 
best fit values of the joint fit between A(S\Lk), A(h\Lk) and 
A(iy|L fc ). 



FIG. 6: Height distribution A(h\Lk) scaling collapse 

(color online) We collapse the simulation data using the uni- 
versal exponents of r = 1.24, Vk — 0.45, and £ = 0.62, the 
best fit values of the joint fit between A(S\Lk), A(h\Lk) and 
A(w\L k ). 



Finally, the respective scaling collapses for the size, 
height, and width distributions are shown in Figures [5j|6j 
and [7| Although scaling collapses are very useful in veri- 
fying critical behavior, we argue that they may be prob- 
lematic for the purpose of determining critical exponents, 
and one should fit and make use of functional forms. In 
Appendix |A 5| we will show how scaling collapses are un- 



able to incorporate the effects of corrections to scaling, 
and how these corrections may cause a drift in the critical 
exponents. 
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parameter 


best fit 


standard errors 


systematic errors 






in linear approx. 





Universal Exponents 



Shared Exponents 



T 




1.2414 


± 0.0006 


± 0.04 






0.4513 


± 0.0001 


± 0.008 


c 




0.6155 


± 0.0004 


± 0.02 


A,Sk(Sk) = 


eMUsSt /2 -Z s S 5 k s ) 






Us 




0.173 


db 0.003 


± 0.2 


Z s 




0.0099 


db 0.0002 


± 0.01 


8s 




1.832 


db 0.004 


± 0.3 


Ahk(hk) = 


eMUhhl /2 - Z h h{ h ) 






u h 




0.94 


± 0.01 


± 0.9 


z h 




0.307 


± 0.004 


± 0.3 


5 h 




1.255 


± 0.003 


± 0.2 


A w k(w k ) = 


- exp(U w wl /2 - Z w w s k w ) 






u w 




0.401 


db 0.005 


± 0.6 


Z w 




0.0291 


db 0.0004 


± 0.2 


S w 




2.202 


db 0.005 


± 0.9 



Non-Universal Exponents 



exp(Af/x + A^/x 2 ) 



At 


-0.36 


± 0.03 


± 2 


A s 2 


-0.35 


± 0.06 


± 4 


A\ 


2.30 


± 0.04 


±2 


A\ 


-1.90 


± 0.04 


± 2 


A? 


-0.99 


± 0.03 


± 2 


A2 


-0.06 


± 0.03 


± 1 



TABLE I: Best Fit exponents and parameters Here are the results of our joint fit for the size A(S\Lk), width A(w\Lk), 
and height A(h\Lk) distributions. The corresponding universal scaling forms which were fit are quoted in the table alongside 
the parameter results; on the bottom of the table are multiplicative corrections for each distribution, with x equal to either 
S, w, and h. Here systematic error bars which account for errors in the theory (see Appendix | A 4| for explanation) are given. 
The traditional standard error bars are typically ~ 64 times smaller than the systematic error bars quoted; however, they are 
a gross underestimate of the actual errors expected since our theory is both highly nonlinear and sloppy 35 . We quote each 
parameter to the significant figure indicated by its standard error, since the parameters are strongly correlated, truncating each 
parameter to its significant figure would yield a poor fit. 
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FIG. 7: Width distribution A(w\Lk) scaling collapse 

(color online) We collapse the width distribution data using 
the universal exponents of r = 1.24, Vk — 0.45, and f = 0.62, 
the best fit values of the joint fit between A(S), A(h) and 
A(w). 



maximum height of an avalanche, but the distribution 
of heights given by random cross sections of avalanches. 
Let A(h x \Lk) dh x be the fraction of the system area con- 
sisting of points (x, y) where the enclosing avalanche 
has vertical cross-sectional height at x in the range 
(h x ,h x + dh x ), then, 

A(h x \L k ) ~ (/iJ4) (2 - t)(1+c) / c A^(^/4)A- (15) 

This distribution gives a different measure of the typical 
shape of an avalanche. The local height distribution is 
shown in figure [8) Here the curves show a fit we have 
generated with a scaling function of the fit to 11 spanning 
avalanches (Eq|32| which cross both window boundaries, 
taking the limit ofW = 0. Our best measure of the local 
heights is equivalent to a 11 distribution with window 
size 1. Details for how this works are provided in the 
section IV CI 



E. Joint distributions and multivariable scaling 
functions 



10" J 



1(T 



IO" 4 



io- 




• • L,k = 16384, IO" 7 
a a L,fc = 16384, 10 6 
▼ ▼ L,fc — 8192, 10~ 5 



< < L,fc=4096, 10 4 
> > L, A; =4096, 10 3 
■ ■ L,k =4096, 10~ 2 



10' 



10 



10 Z 



10 J 



FIG. 8: Area-weighted local height distributions (color 
online) Here are the area- weighted local height distributions, 
the fraction of area taken up by a cross sectional height h x . 
The fits shown in the figure were with the form of Eq. |15| 
where Ah x k is the scaling function of the fit to 11 spanning 
avalanches (Eq: 32 ) which cross both window boundaries, tak- 
ing W = 0. Details for this function are explained in sec- 
tion M 



D. Local height distributions 

In our analysis of the next section, we will make use 
of another scaling function of the same form: not the 



Once we have distributions for the measures of size 5, 
width w and height h, we can also explore the forms 
of joint distributions. The area A(w, S\Lk) dS dw of 
avalanches in the range of size S + dS) and widths 
(w,w + dw) will go to zero strongly if the size S becomes 
either much larger than or much smaller than the typical 
size w 1 ^^ of an avalanche of width w - so we may fac- 
tor out any combination of powers of w and S without 
changing the singularity. It still makes sense, though, to 
factor out the Lfc-dependence. If we choose to factor out 
powers of S, therefore, we find 



A(w,S\L k ) 
( S \ 



(16) 



L 



i+C 



A 



wSk 



w w 



L (r-2)(1+C) ^l-r-l/Cl+C)^ 



wSk 



w 



w 



where again we have a power of a dimensionless scal- 
ing variable, times a scaling function, divided by Sw ~ 
551/(1+0 s i nce A(w,S\Lk) is multiplied by dS dw in its 
invariant form. 

The last joint distribution that will be useful is re- 
lated to right-most pieces of an avalanche. Consider the 
right-most piece of width x of an avalanche of total width 
w and total size S; let this segment have size s. (This 
will be the size measured by a window that cuts the 
avalanche at the left-hand window boundary at x.) Let 
A(s,w,S\x,L) be the fraction of the system covered by 
such avalanche pieces. Then, in the same logic as before, 
this five-variable distribution can be written as a power 
law times a universal function of four variables: 
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A(s,w,S\x,L) L 



(r-2)(l+C) c 



1/(1+0^4 



swSxk 



s V(i+C) ' s ' #V(i+C) ' L fc 



(17) 



r 



One can clearly work out scaling forms for joint distri- 
butions of several variables and other combinations. The 
ones we have discussed here will be needed in our analysis 
of windowing effects. 



V. WINDOW EFFECTS 

Now that we have laid the groundwork for explor- 
ing the shapes of avalanches, we focus on analyzing 
avalanches inside a viewing window. In this section, we 
focus on how to define the right power laws and scaling; 
we also give results for fits, extracting critical exponents. 
In the next section, we go into more detail about the scal- 
ing shapes of these distributions - the universal scaling 
functions for avalanche sizes viewed through windows. 

In imaging experiments one often runs into the prob- 
lem of not being able to see the whole system, dis- 
torting the avalanche size distribution. In particular, 
for Barkhausen noise, typical magnetic avalanches span 
many decades in size, far beyond the spatial resolution of 
optical microscopes. The natural solution is to take mea- 
surements at a variety of magnifications. Even though at 
the weakest magnifications the window size W > and 
most avalanches avoid the window boundaries, the ef- 
fects of the boundaries will always dominate at the high- 
est magnifications. The analysis in this section not only 
provides a method to correct for finite-size-like window 
effects on exponents, but allows us to actively make use 
of all the data for a range of magnifications. 

We show in detail in this section and the next how 
the avalanches which cross different boundaries exhibit 
distinctly different size distributions and critical expo- 
nents (Figure [lj. As described in section \H\ we consider 
the avalanches measured in an infinite strip of width W 
(Figure [TJb)), for a system with characteristic length L^. 
We separate avalanches into different categories: inter- 
nal avalanches (00), split avalanches (10 and 01), and 
spanning avalanches (11). Let us call Aoo(s| W, Lk) ds the 
area fraction covered by such avalanches with sizes in the 
range s, s + ds. (For Aoo, the segment size equals the to- 
tal size.) The split (10 or 01) avalanches will have area 
fraction Aio(s|W, Lk) ds for each s. The distribution Aqi 
(Figure [TJf)) of avalanches touching the right boundary 
naturally equals Aio on average. Au(s\W, Lk) ds is the 
fraction of the strip spanned by 11 spanning avalanches. 
We mentioned in section [ill that the 00 avalanches have 
a power law that matches the full system, whereas the 
10, 01, 11 avalanches all have modified power laws with 
a smaller exponent r. 

Besides different power law scaling, the universal scal- 



ing functions for these different avalanche distributions 
are also distinct. In particular, the cutoff dependence on 
window size is different for internal avalanches and split 
avalanches, while the spanning avalanches have both an 
outer and an inner cutoff due to the window size (since 
avalanches must be large enough to span the window). 
We present the fits of these universal scaling functions in 
this section and discuss their shapes in more detail in the 
next. 

We know that all avalanches in the window are of one 
of the 00, 10, 01, 11 types, so 



/ 



dsA 00 (s\ W, L fe ) + 2A w (s\W, L k ) + A u (s\W, L k ) = 1. 

(18) 

As in the previous section, consider how each of these 
distributions A zz rescales under a coarse-graining by a 
factor b. The zz denote our indices for the various win- 
dowed distributions (00, 10, 01, 11). Each A zz ds, being a 
geometrical quantity (a fractional area), must be invari- 
ant under rescaling (with two invariant scaling variables) : 

A zz (s\ W, L k ) ds = A zz (s/b 1+ <\W/b, L k /b)ds/b 1+ < 

= (l/s)B zz (s/L k +c ,W/L k )ds. (19) 

However, this is clearly not the form which makes the 
size and window-width dependence of the avalanche sizes 
manifest. We are allowed to factor powers of the invariant 
scaling variables W/Lk and s/L 1 ^^ out of the scaling 
function B: 



A zz (s\W k ,L k ) 



1 



s \ L 



i+C 



s W 



-(T«-2)(l+C)+*>**TI/-t>« £ .l 



W- V **8 l - T "A Z 



W 



(20) 



The appropriate powers of s, and W to pull out- 
side depend upon which of the three distributions we are 
considering. For the distributions Aqq and Aio, ^oi? we 
choose t zz and v zz (powers of the invariant scaling vari- 
ables X = s/L 1 ^^ and Y = W/Lk) to make the result- 
ing scaling function go to a constant at small X and/or 
Y. This way the power laws we pull out describe the 
behavior of the limit of small s, and the way in which 
the avalanches are cut off by the window size (as s ap- 
proaches W 1+ ^) are described by the scaling function. 
On the other hand for An there are no small avalanches 
(they have to be large enough to span the window), and 
as Lk — > oo, all avalanches span the window and become 
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11 avalanches, and in this limit the distribution will not 
go to zero, so for this distribution we instead pull out 
powers of W and s. 



touch the boundaries. This can be computed explicitly 
from the function A(w,S\Lk) (Eq. 16) which gives the 
area covered by avalanches of width w and size S (note 
that for 00 avalanches, the segment pieces s = S): 



Internal Avalanches 



First let us consider Aqo(s\W, Lk)ds, the window area 
spanned by avalanches of sizes in (5,8 + ds) that do not 



A o(s\W,L k ) 



w 



I dw 

J a 



w 



w 



W -w 
W 

>_ l (t-2)(i+o q i 



A(S,w\L k )dw 

51— mi+o AwSk 



-(r-2)(l+C)ol 

J I. & 



_1/(1+C) [ dW ^ A ^ 



w 



w 



5^(1+0 ' L k 
w w 



(21) 



r 



where (W — w)/W is the probability that an avalanche 
whose center lies within the window is entirely contained 
in the window (i.e., the avalanche center lies within (W — 
w)/2 of the center of the window). Changing variables 
from w to Q = w/s 1 ^ 1 ^ , 



r W/s~' - 

A 00 (s\W,L k ) = L^- m+O s 1 - T / 
s i/(i+0« 



1/(1+0 



dn 



(i + C) 



W 



^(r-2)(l+C) 1 



rl+C 



5 1 -M 00 (5/W 1+c ,Ty/L fc )(22) 

with no explicit dependence on the window width W 
(so voo = 0), and the same critical exponent Too = r 
that is given by the non-windowed distribution. Note 
that Eq. ([22]) is of the general form given by Eq. (19) 



and ( 20 ). This scaling equation is also consistent with our 



numerics: the normalization of the distribution for small 
avalanches is independent of W, and Too = 1.26 ± 0.02 is 
consistent with the bulk r = 1.24 ± 0.04. 

Using this scaling form, we fit the 00 data (jointly with 
11 and 10 data) with a scaling function given by a pa- 
rameterized functional form: 



A 00 (s\W,L k )=L 
Zoos s k 00 +C o( 



(r-2)(l+C) 1_ T , ,rp 
Sk 



TT 1/2 



W noo 



) m °°)eM(Ar/s + A™/s 2 )) 



(23) 



with Sfc = s/L^l + C)- Out fit for the parameter noo is 
1.62, which we believe to be 1 + C (see next section). This 
makes the term s k /W£ 00 = S/W 1+( * which is another 
natural invariant scaling variable. 



Figure [9] shows the results of a nonlinear least squares 
fit, with shaded areas as estimations of fluctuations in 
the theory corresponding to systematic error bars on our 

C, (and v k which is hidden 



In Eq. 23 



parameters. 

in L k ) are universal exponents shared amongst the three 
different distributions, A® and A^ are (non-universal) 
analytic corrections to scaling reflecting lattice effects on 
small avalanches, and the other parameters encapsulate 
the shape of the universal scaling function *4 o- The fit- 
ted results for the other universal and non-universal pa- 
rameters are quoted in Table [TTJ We describe the scaling 



shapes and their motivation in more detail in section VI 



One may note the exponent £ in our fits is fixed to the 
literature value of 0.63. If we allow for a free fit on all 
the parameters, it shifts to ( = 0.68 ±0.02. Although the 
free fit is 2.5a away from the accepted value of £, the fit 
with fixed £ = 0.63 has only a 50% higher cost than the 
free fit minimum, suggesting an average of 0.5<j drift on 
the parameters instead of 2.5a as seen in Table \U\ This 
suggests three cautions (1) The estimate on our system- 
atic error (0.02) is a lower bound estimate, and in fact 
the systematic error should be higher. (2) The fact that 
our systematic error should be higher also implies that 
the scaling functions are imperfect and may be improved 
upon. (3) There could be corrections due to a crossover 
that depends on both A and k which we have not ac- 
counted for, which are distorting the fit. The subtleties 
and nuances in measuring the exponent and the pos- 
sible origins of this drift are discussed in more detail in 
Appendix [Bl 
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B. Split Avalanches 



Next, consider the avalanches that are split by one 
side of the window, say the left side, with the distri- 
bution Aio(s\W,Lk). Physically, for small avalanches s 
and large L k /W this is clearly proportional to 1/W: the 
small avalanches extend only a small distance into the 
window (smaller than the window width), so the frac- 
tional area covered by them is proportional to one over 
the width of the window. This leaves us with a scaling 
form 





= 10 


5 ,0.36 


< < k,W s 


= 10 


6 ,1.02 


o o k,W s 


= 10 
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A w (s\W,L) = 

J_ r (r'-2)(l+0 ,l-r' + l/(l+0 A 

w s A 



(24) 



10 



(s/Ll + t,W/L k ). 



FIG. 9: Internal Avalanches Data and Fit (color online) 
Shown here are the area-weighted size distributions for inter- 
nal (00) avalanches. The lines are the joint best fit of Aoo, 
Aio and An to the functional forms of equations |23| [26] and 
[32] whereas the shaded areas are the fluctuations in theory 
corresponding to the systematic covariant errors on our ex- 
ponents and parameters (individual parameter best fit values 
and errors are quoted in Table |lT| . 



with r' to be determined. Note again that Eq. (24) is of 
form Eq. (19). 



We can also write A\§ in terms of the distribution of 
right-most pieces A(s,w, S\x, L k ) from Eq. [Tt| integrat- 
ing over all possible sizes 5, all possible widths w, and 
all possible pieces x (from lattice size a to window width 
W): 



f°° f°° f w dx 
A 10 (s\W,L k )= / dS dw —A(s,w,S\x,L k ] 

J a 2 J a J a 

( s Y~ T i r Ll * +< ,„ r L - , r w dx 



2 ol/(l+C) 



s 



1 



"4 s - 



dS I dw 

J a 



As 



w ^swSxk I i/(i +C ) > s > 51/(1+0 ' L k 



d ( w 

a/sV(i+C) VsV(l- 



+0 



/>w/ s 1/(1+c) s i/(i+0 / x x / x 

x y a /,i/(i + o ~^^ d C^vaTcy; AswSxk 

= ^4 T - 2)(1+C) ^- r+1/(1+C) Ao(,/4 +C , W/L k ) 



w 



w 



+0' s '5V(i+C)'L fc 



(25) 



cte/W is the relative probability that the avalanche inter- 
sects the left-hand boundary, and we have changed the 
integration limits at 00 to the avalanche length scale L k . 
(For w < x < W, the original distribution is naturally 
zero.) After we rewrite the integration variables in terms 
of the invariant scaling variables, we can organize the 
form of the scaling function into the form of Eq. [24j This 
tells us that r' = r. These results are consistent with our 
numerical fits: W has an exponent of minus one, and r' 
is equal to the system r. 

With the correct power laws pulled out, we can now 



write down a function to describe the data and cutoff: 
A 10 (s\W,L k ) = ^Lt 2)(1+ °s 1 -^ 1 ^ 



cxp(-(T 10 + U w s\r + Z 10 s s k " + C 10 (^) mi °) 



4 /2 



eM(A ( r/s + A^>/s')) 



(26) 



Again, our best fit nio is nearly 1 + £, so Sk/W^ 10 
s/W 1 ^^. Also note that this distribution has the same 
functional form as Aqq in eqn 



,i/(i+C) 



/W in front. 



23 aside from a factor of 
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Figure [10] shows the results of a joint nonlinear least 
squares fit with the 00 and 11 avalanche data, with 
shaded areas representing estimations of fluctuations in 
the theory corresponding to systematic error bars on our 
parameters. Here, as in the 00 distributions, r, £ (and v k 
which is included in Lk and the scaling variables Sk and 
Wk) are universal exponents shared amongst the three 
different distributions, A\° and A^ are (non-universal) 
analytic corrections to scaling reflecting lattice effects on 
small avalanches, while the other parameters describe the 
shape of the universal scaling function Aio(sk, Wk)- Fit- 
ted results for the other universal and non-universal pa- 
rameters are quoted in Table [TTJ We describe the scaling 
shape and its motivation in more detail in section VI 



the avalanches have size s — Wh x , where h x is the height 
of the vertical cross section of the avalanche. Hence we 
can write An in the limit W — >• in terms of the distribu- 
tion A(h x \Lk) of randomly chosen vertical cross sections 
of avalanches (eqn [l5|) , choosing h x = s/W: 

Au(s\W=0,L k )ds ~ w ^ dh x A(h x =s/W\L k ) 

A 11 (s\W=0,L k ) ~ w _> ±A(h x =s/W\L k ) (27) 



where dh x /ds = 1/W. Remembering from Eq. ( [15] ) 

that A(h x \L k ) ~ (^/4) (2_r)(1+C)/C A cc fe(^/4)/^ 5 
we substitute s/W for h x and take the limit of Wk — > 
in Eq. [27] to give: 




io u io 1 i(r io j io 4 io' io b 



Au(s\W = 0,L k ) ~w k ^o 

\ (2-r)(l+C)/C 



1 



w ywLi, 



Ah x k 



/(s/w)m 



we cancel the two W's and add the dependence on the 
second scaling variable Wk = W/Lk to derive the scaling 
form for An: 



A u (s\W,L k ) 



(2-r)(l+C)/C 



WL\ 



An (s/(WL c k ),W/L k y (29) 
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Here limy^o An(X, Y) = Ah x k(X) and thus 
J Au(X, 0) dX = 1 (implied by the fact that al- 
most all avalanche area touches both boundaries as 
W/Lk —> 0). Also notice that since: 



FIG. 10: Split Avalanches Data and Fit. (color online) 
Shown here are the area-weighted size distributions for split 
(10) avalanches with different k and window size W. The 
lines are the joint best fit of Aoo, Aio and An to equations 
[23] [26] and [32] whereas the shaded areas are the fluctuations 
in theory corresponding to the covariant systematic errors on 
our exponents and parameters (individual parameter best fit 
values and errors Table [5}. 



v s s L k 1 ( v 

X = WZi = I^ow =Sk -W m 

we can rewrite Au(s\W, Lk) as: 

1 / s t \ (2 " r)(1+c)/c ( Sk 



(31) 



C. Spanning Avalanches 

Finally, consider the spanning avalanche distribution 
An. First, remember that most of the area in gen- 
eral is spanned by the largest avalanches (since r < 2). 
Therefore, as Lk — » oo, 100% of the area is covered by 
avalanches of widths much larger than W and hence An 
must integrate to one in this limit. This makes it natural 
to pull out only powers of W and s outside the scaling 
form for An. 

Second, notice that the size of an 11 avalanche is ba- 
sically given by its height. More specifically, as W — » 0, 



Figure [TT] shows the results of a joint fit with the simu- 
lation data of the previous 00 and 10 distributions. For 
the 11 distributions we use the functional form: 



A n (s\WL k ) 



1 



Sk \ 



(2-r)(l+C)/C 



exp(-(T n + Uus 1 ,/ 2 - 



Znst+Dn 
exp(A^ 11) /s) 



Sk 



Sk 



-ml 



(32) 



r (and also ( and v k which are hidden in the scaling 
variables s k and Wk) are universal exponents shared 
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amongst the three different distributions, A\ x is the (non- 
universal) analytic correction to scaling reflecting lat- 
tice effects on small avalanches, while the other parame- 
ters describe the shape of the universal scaling function 
Wfc). Note that we don't include the term A\ Xl 
as we have done in the 00 and 10 distributions; this term 
turns out to be the same as another term in the univer- 
sal scaling function in the l imit of W — >• 0, and so it is 
redundant. (See Appendix A 5). The best fit universal 



and non-universal parameters are given in Table [II] We 
discuss in more detail the motivation and form of the 



scaling function An(X,Y) in section VI 
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FIG. 11: Spanning Avalanches Data and Fit (color on- 
line) Shown here are the area-weighted size distributions for 
spanning (11) avalanches with different k and window size 
W. The lines are the joint best fit value using the functional 
forms using equations [23] [26] and [32] whereas the shaded 
areas are the fluctuations in theory corresponding to the sys- 
tematic errors on our exponents and parameters (individual 
parameter best fit values and errors are quoted in Table [lT| . 



We also test the limiting case of our scaling function 
with our data sets of W = 1 in Figure [8] of section |IVE 



The curves drawn in Figure [8] are with the function given 
in equation [32] using the best fit values of the joint fit of 
An, A10 and Aqq. Notice that the predictions of equa- 
tion [32] matches the data for Wk = 0, indicating our 
function satisfies limy Ai i(X, Y) 



A hxk (X) as ex- 



pected. 



VI. SCALING SHAPES AND RESULTS 

In the previous section we wrote down scaling functions 
for each type of avalanche inside a window. In principle 
there are many possible parameterizations we can choose 
that would be able to capture the behavior of the data. 
In this section we explain how and why we chose each 
one, and also discuss the scaling function in the limit of 
small windows. 



A. Scaling shapes and functional forms 

We would like to capture the scaling behavior of both 
the finite size of the avalanches, and the effect of the 
window size on the distributions. We choose forms in- 
spired by a functional renormalization group expansion 
for static avalanche size distributions for all universality 
classes [31] , and further motivated by heuristic arguments 
for the cutoff dependence on the two scaling variables, 
s k = s/L 1 ^ and W k = W/L k . 

Similar to in section [IV C[ we start with an avalanche 
size distribution of the form: 



We expect that as W — » 00, we will not see the effect of 
the window size, and the scaling forms will go to the limit 
of our proposed avalanche size distribution in Eq. 33 
Keeping this in mind, we would like to write a function 
of A(sk, Wk), including the effect of the window size. For 
the 00 and 10 distributions, we expect the avalanche sizes 
to be cutoff by the window size W when W/Lk < 1, hence 
a cutoff dependent on Sk/Wg should be expected, where 
n > 0. While when W » Lk, their scaling forms should 
go to the limit of our proposed avalanche size distribution 
in Eq. [33] Therefore for the 00 and 10 avalanches, we 



propose a scaling form of: 

A zz (s k ,W k ) 

exp(-(T zz + U zz sl /2 + Z zz s* k * 



■c z 



Sk 



WV 



A(s k ) = exp(-(T + Us]! 2 + Zs{)) 



(33) 



We have heuristic arguments for what this n zz value 
should take. Figure [12] a schematic of an avalanche cut 
by different windows, is meant to help illuminate our 
discussion. For the 00 internal avalanches, the largest 
avalanche contained fully within the window should have 
a width w that is roughly W. And since s max ~ 

Wmaxhmax ~ ™lrtaii then S max ~ W 1+C and it fol- 
lows that Sk ~ W£ + ^ , giving us an expectation value 
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parameter 


best fit 


standard errors 


systematic 


drift from best fit 






in linear approx. 


errors 


with free f 


Universal Exponents 


T 


1.2636 


± 0.0006 


± 0.02 


-0.006 




0.4630 


± 0.0002 


± 0.01 


-0.02 


c 


0.63 (fixed) 


± 0.0007 


± 0.02 


+0.05 


Aoo(s k ,Wk) 


exp(-(T 00 + Uoos\ /2 + Zoos{ OQ 


I fit s k \m 00 \ 
k 


Too 


2.488 


± 0.004 


± 0.1 


-0.01 


Uoo 


-0.150 


± 0.005 


± 0.1 


+ 0.04 




0.0040 


± 0.0004 


± 0.01 


-0.0009 


Soo 


2.21 


± 0.03 


± 0.9 


+0.06 


Coo 


5.60 


± 0.01 


± 0.7 


+1.8 


moo 


1.371 


± 0.003 


± 0.1 


-0.04 


noo (1 + C) 


1.621 


± 0.004 


± 0.7 


+0.04 


Aio(s k ,Wk) 


exp( — (Tio + Uios\ /2 + ^ios^ 10 


+ ^10(^0 ) mi °) 

k 


Tio 


1.437 


± 0.004 


± 0.1 


-0.1 


Uio 


0.244 


± 0.244 


± 0.1 


-0.03 


Z\o 


0.027 


± 0.001 


± 0.03 


+0.005 


$10 


1.64 


± 0.01 


± 0.4 


-0.06 


Cio 


1.153 


± 0.004 


± 0.2 


+0.7 


mio 


1.962 


± 0.005 


± 0.2 


+0.04 


nio (1 + C) 


1.624 


± 0.004 


± 0.1 


+0.06 



Au(s k ,W k ) = 

exp(-(Tii+f/ii4 /2 + Zii4 11 



+ D 11 (^r 1 



_ ^ -m2j 





0.47 


± 0.03 


± 1.2 


-0.3 


t/n 


-0.5 


± 0.1 


± 3.6 


-0.5 




0.21 


± 0.06 


± 1.7 


+0.4 


*ii 


1.102 


± 0.03 


± 1.0 


-0.12 


On 


0.52 


± 0.03 


± 1.0 


+0.1 


Cn 


0.83 


± 0.05 


± 1.6 


-0.3 


ml 


1.48 


±0.01 


± 0.4 


-0.0008 


m2 


1.64 


± 0.02 


± 0.6 


-0.02 


nii(l + C) 


1.655 


± 0.004 


± 0.1 


+0.02 


Corrections to Scaling Af z / 


s + A z 2 z /s 2 






A , 


-0.94 


± 0.02 


± 0.5 


+0.06 


Af 


0.27 


± 0.01 


± 0.4 


-0.04 


A\° 


0.15 


± 0.01 


± 0.4 


-0.2 


Af 


-0.07 


± 0.01 


± 0.3 


+0.1 


Al 1 


0.8 


± 0.07 


± 2.1 


-0.3 



TABLE II: Best fit exponents and parameters for windowed distributions. Here are the results of our joint fit for 
the windowed Aoo, An, A±o distributions. The corresponding scaling forms which were fit are quoted i n the table alongside 
the parameter results. Here systematic error bars which account for errors in the theory (see Appendix A 4 for explanation) 
are given. The traditional standard error bars are typically ~ 30 times smaller than the systematic error bars quoted. The 
last column is the drift in parameters seen when allowing ( to be a free parameter. Notice that these numbers are more 
or less similar or smaller than the estimated systematic error, except for £. (The problems in measuring £ are discussed in 
Appendix [bJ As in Table [I] we quote several digits more than the error bars warrant for individual parameters, because the 
errors are strongly correlated; truncating each parameter to its significant figures would yield a poor fit. 
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FIG. 12: An avalanche cut by windows in the extreme 
limits, (color online) Drawn here are two fronts separated 
by an avalanche event. Here we are depicting cases where 
this avalanche (or its segment) is the maximum avalanche 
size for the 00, 10, and 11 at different window widths. Boxes 
of different widths and colors are used to show the cases in 
which this may happen. The main avalanche may count as 
a 00 avalanche for a wide window, while part of it would 
count as a 11 avalanche for a smaller window; it could also 
count as a 10 avalanche if it happens to cross the window 
boundary. This figure illustrates our arguments for the shape 
of the cutoff (the exponent n zz ) given a window size W for 
the 00, 11, 10 cases. Another small avalanche is drawn for 
the 11 case to show that the minimum size to cross the 11 
window also introduces a separate cutoff. 



of rioo = 1 + C Numerically, we find noo ~ 1-62 when we 
fix ( = 0.63. The fit plotted against one of the scaling 
variables and also the contour plot of the scaling shape 



is shown in Figure 13 



For the 10 or 01 spilt avalanches, since we are effec- 
tively measuring the ends of avalanches that spill into 
the window, nio depends on what the shape of the 
avalanche is at the edges. The largest portion of an 
avalanche to spill into the window will again be limited 
by the size of the window W. Here the shape follows 
the roughness of the two fronts preceding and following 
the avalanche, where h(x) ~ for each, so plausibly 
hedge = hafter — hb e f ore ~ W** . The size is then limited 
by s max ~ w max h edge ~ W 1+c , giving us an expectation 
value of nio = 1.63. Numerically, we find nio ~ 1-62 (( 
is estimated from 0.62 — 0.72 in our various measures), 
matching our expectation. The fit plotted against one of 
the scaling variables Sk and also the contour plot of the 



scaling shape is shown in Figure 14 

Now we move on to discuss the 11 spanning avalanches. 
Here the situation is slightly more complicated than the 
previous two cases, due to the distribution being strongly 
cut off at two length scales, as one may note from the 
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FIG. 13: Internal Avalanches Scaling Function, (color 
online) (a) Scaling collapse showing Aoo(sk, Wk) as a function 
of Sk- The separate curves show the effects of the scaled 
window size Wk. (b) Logarithmic contour plot of best fit 
scaling function against both scaling variables Sk and Wk- 
Each contour reflects a drop of a factor of e in the scaling 
function. The black dots are at locations of the simulated 
data points used in the fit indicating where the fit should 
be a reliable prediction. The red solid line is log 10 Wk = 
noo log 10 Sk which is the slope at the large avalanche cutoff, 
with noo = 1.62, the best fit value, (noo = 1 + C — 1-63 is the 
expected value from our heuristic arguments). 



shapes of the distributions shown in Figure 11 First of all 
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FIG. 14: Split Avalanches Scaling Function, (color on- 
line) (a) Scaling collapse showing Aio(sk,Wk) as a function 
of Sfc. The separation between curves shows the dependence 
on the scaled window size Wk- (b) Logarithmic contour plot 
of best fit scaling function against both scaling variables Sk 
and Wk- Each contour reflects a drop of a factor of e in the 
value of the scaling function. The black dots are at locations 
of the simulated data points used in the fit. The red solid 
line is log 10 Wk = nio log 10 Sfc which is the slope at the large 
avalanche cutoff, with mo = 1.62, the best fit value. 



the avalanches need to be large enough to cross the win- 
dow, implying an inner cutoff that depends on Wl jsk 



WjJ +< */sfc); here the argument for the minimum size 11 
avalanches follows from a similar argument for the max- 
imum size 00 avalanches. The smallest avalanche that is 
able to span the window will have a width w m i n = W . 
whereas h min ~ w£- , and s min ~ h m i n w m i n ~ W 1+c > . 



However, in this case s k should be in the denominator 
of the scaling variable, since for s max » s > s m i n the 
probability of the having a spanning avalanche grows as 
s increases. For the outer cutoff we expect that the max- 
imum size is given by the window size W and the typi- 
cal maximum height, i.e. s max « Wh rnax which implies 
that s k max « W k max (h max /L c k ). This rescaled height 
(h max /Ly is constant. Therefore the cutoff for the 
large size avalanches should depend on Sk/Wk- Hence, 
we propose the scaling form below for the 11 spanning 
avalanches: 



An(s k ,W k ) = exp(-{T n ■ 

-ra 2 

/ .s ?, \ 



Sk 



U n s k 



1/2 



Sk 



)) (35) 



where Cu controls the strength of the inner cutoff and 
Dn the outer cutoff. 

Figure [l5^a) shows the shape of the scaling function 

lb) 
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(i.e. the cutoff is for s m i n /s < 1 and so 



Js 



plotted against one scaling variable and Figure 
gives the contour plot of this function against both vari- 
ables. The best fit value of n\\ is n\\ ~ 1.65, whereas 
the expected was n\\ = 1 + ( = 1.63. 

Finally, one may note that in our system all the n zz 
turn out to fit to our expected 1 + ( within the error 
bars of £. One may be tempted to set n zz as 1 + £ and 
have fewer fit parameters in one's form; however, we rec- 
ognize that our geometrical arguments do not hold for 
front propagation that is super-rough with £ > 1, or for 
models that allow overhangs, so we leave ti zz as a free 
parameter to signify this geometrical constraint. 



B. The limit of small windows 

Although there are noticeable imperfections in the the- 
ory function, the agreement is impressive between theory 
and simulation as seen in figures [l3j [l4j and [15] . 

The scaling function for each of the three distribu- 
tions is a competition between two types of terms: the 
rescaled size s k = S/L 1 ^^ and the rescaled window size 
Wk = W/Lk. Upon examining the fit parameters of the 
scaling function, all three distributions share the charac- 
teristic that at small W k the terms with pure powers of 
s l J 2 and s 6 k zz become unimportant, leaving only s k /W k zz 
for the 00 and 10 distributions, and for the 11 a s k /Wk 
term. Notice that for the 00 and 10 distributions, since 
n zz is 1 + the L k dependence disappears for the univer- 
sal scaling function at small Wk- Therefore, the shape of 
the scaling function is cut off mainly by the window size 
W. In fact removing the s k terms for these functions does 
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FIG. 15: Spanning Avalanches Scaling Function, (color 
online) (a) Scaling collapse showing An(sk, Wk) as a function 
of Sfe. The curves move leftward and become more sharply 
rounded with increasing Wk- (b) Logarithmic contour plot 
of best fit scaling function. Each contour reflects a drop of 
a factor of e in the scaling function. The black dots are at 
locations of the simulated data points used in the fit. The 
upper orange solid line is the slope of the contour plot at the 
small avalanche cutoff, and has log 10 Wk = 1/1.65 log 10 Sk 
where 1.65 is the best fit mi value. The orange dashed has 
\ogWk = 1/1.63 log 10 Sk where 1.63 is the 1 + C value. The 
lower red solid line is log 10 Wk = log 10 Sk which is the slope 
at the large avalanche cutoff. 

not affect the shape for Wk < 5. For experiments that 



study systems in the same universality class as this one, 
this implies that data may be measured at large magnifi- 
cations (small windows) and fit to extract exponents and 
scaling behavior without the extra, often unknown, scale 
of L k . 



VII. SUGGESTIONS FOR EXPERIMENTS AND 
CONCLUSIONS 

What does our analysis imply for current experiments? 
How should one conduct the experiment and analyze the 
data? Here we discuss for the particular case of magnetic 
avalanches in Barkhausen noise, how to take into account 
window effects and further enhance the collection of data. 

There are two optical methods for detecting avalanche 
distributions for Barkhausen noise in 2D thin films, 
and both make use of the magneto-optical Kerr effect 
(MOKE). When a polarized beam of light reflects off a 
magnetized sample, the reflected polarization is affected 
depending on the magnetization. A second polarizer can 
be used to filter this signal, and then using either a pho- 
todiode [16] , or an optical microscope |5[ |6] , we can col- 
lect data about the avalanches from the signal. For ex- 
periments using a photodiode (let us call this "laser re- 
flectometry" ) , one can only measure the total magneti- 
zation change over time, and not individual avalanches. 
For experiments using an optical microscope (let us call 
this "avalanche visualization"), one can resolve individ- 
ual avalanches and their shape. 

In laser refiectometry experiments [16] we only have in- 
formation for the magnitude of magnetization as a func- 
tion of time, and cannot see which avalanches touch the 
boundary. Furthermore, in current techniques the laser 
spots are Gaussian in shape, and do not have sharp 
boundaries. However, there seems to be no fundamental 
reason why the illuminated region could not be optically 
generated with uniform illumination and sharp edges, up 
to some diffraction limit depending upon the geometry 
of the experiment. (A typical avalanche of interest is a 
few microns in size, large compared to the wavelength of 
optical light which is 400-700 nm). If one could make 
the edges of the laser spots sharper, one could adjust 
the laser spots to flicker between two sizes, one with a 
radius slightly larger than the other. Events that occur 
with the same magnitude in both the large size mea- 
surement and the smaller size measurement would be 00 
internal avalanches. More elaborate sequences of spot 
shapes could be used to further distinguish 01 and 10 
split avalanches from 11 spanning avalanches. 

For avalanche visualization experiments [5j |6j [22] , it is 
straightforward to separate the data into 00, 10, 01, and 
11 avalanches for systems which have our strip geometry 
and fiat fronts with £ < 1. Our analysis will remain valid 
with minor corrections due to real-world experimental 
circumstances. For example, sometimes the propagation 
direction of the front is not parallel to the top and bot- 
tom boundaries. In this case there would be corrections 
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depending on the angle of the tilt and the size of the 
window. The local heights h x would need to be adjusted 
with the factor of cos 0. The 10, and 01 avalanches would 
be cut off at an angle, but for self-affine (short, wide) 
avalanches these size corrections are likely irrelevant. 

For materials with dipolar interactions and zig-zag 
shaped fronts [5j [36] , we are less confident that our meth- 
ods can be applied without modification. The large ver- 
tical extent of the zig-zag front suggests that all four 
boundaries of the window will matter; therefore we will 
need to divide avalanches into more categories (0000, 
1000, 0100, 0010, 1010...). The analysis for these types of 
avalanches will be more complicated. This would be an 
interesting problem to pursue by simultaneous analysis 
of simulations and the experimental data. 

Knowing how to deal with window effects can be an 
important tool for these visualization experiments. By 
combining data at several magnifications (corresponding 
to different window sizes), we resolve a larger range of 
length scales. Higher magnifications will show the small 
avalanches, while lower magnifications will allow us to 
both capture larger avalanches and explore more fully the 
00 internal avalanches. For example, if our CCD camera 
recording the images has a resolution of 1000 2 pixels, and 
we have a magnification of up to 5x-50x, we can simulta- 
neously explore a range of window sizes and extend our 
effective spatial resolution from 1000 2 to 10000 2 . 

More generally, this paper has provided the tools 
needed to extract from the experimental data for sys- 
tems of similar interface dynamics the critical exponents 
r, £ and the universal scaling functions Aski Ahh, etc. 
For these experiments, we can also measure widths and 
heights and the average shape, giving us an independent 
measure of £, and the universal scaling function Ah x k- 
Our detailed analysis and comprehensive methods pre- 
sented in this paper enable a more powerful interpreta- 
tion of current experiments, and improved construction 
of future experiments. 



Appendix A: Methods and software 

1. Numerical simulation 

The avalanche simulations in this paper were produced 
using a quenched KPZ model [lOj [11], with dynamics 
given by: 



dh(x, t) 
dt 



F-^>+ 7 V 2 /i + A(V/i) 2 +^,/i) (Al) 



where h(x,i) is the height of the front, F the driving 
force, k the "demagnetization field", linear and non- 
linear terms for the KPZ model controlled by the param- 
eters 7 and A respectively, and n gaussian random noise. 
This was run for system sizes of width L 4096, 8192, and 
16384. The simulations have been run in a strip geometry 
(4096 x 8192, 8192 x 16384, and 16384 x 32768) and the 
bottom half of the simulations have been truncated to 



avoid transient effects due to the initially flat front. The 
left-right boundaries have periodic boundary conditions. 

The simulations are done using a discrete cellular au- 
tomaton model, in which the displacement of the string 
h, the time t and the space x are all discretized and 
take integer values [13] [37. For a configuration {hi}, we 
compute the local force Fi at each site z, leading to a 
discretized version of Eq. |A1| 



F i = F-k{h i ) + ^YS- h * 



+nn hi) 
2 



■^T(hi+i - hi^) 2 +r]i(hi) 



(A2) 



where the sum runs over all the nearest neighbors nn for 
the site z, a is the discretization length that we set to 1, 
and rji(hi) is a random force. The automaton dynamics 
are as follows: (1) increment the external field until one 
site is unstable (Fi > 0); (2) determine for each site along 
the interface whether it is stable (Fi < 0) or unstable 
(F{ > 0); (3) advance all unstable sites by one step hi = 
hi + a = hi + 1 in parallel, generating a new value of the 
pinning force rn(hi)\ (4) repeat until no sites are unstable 
(the end of the avalanche); (5) repeat (1-4) until the front 
passes the top of the simulated window. 



2. Nonlinear least squares fitting 

We use nonlinear least squares methods for fitting data 
to theory functions, minimizing a cost defined as: 



C(6) 



1. 



y. 



theory 



y 



data 



(A3) 



Here are the parameters, y the function value, and cr^ 
the error on the data points. The weight a in our case 
is determined by fluctuations from run to run. Namely, 
we bin the data (in equispaced log bins) for each of N 
simulations, calculate the standard deviation of this value 
across runs, divide by ^/N — 1 to get the fluctuation in 
the mean of that bin. [44] 

In our distributions, small avalanches occur more of- 
ten (leading to small error bars), but large avalanches are 
more important, and the smallest avalanches suffer from 
non- universal lattice effects. So during the fitting pro- 
cess there is a tradeoff between fitting the region where 
there is good data and where the variations are most 
important. We use a number of methods to compen- 
sate for this imbalance. (1) We set a minimum error 
bar (1% of the data value) on the data points, making 
the error bars on smaller avalanches larger, and therefore 
decreasing their weight. [45] (2) Analytic and singular 
corrections to scaling can also account for non-universal 
effects. We include analytic corrections to scaling for the 
lattice effects in our scaling functions. These corrections 
appear in all distributions we discuss. A more detailed 
discussion is in section E of this appendix. (3) One may 
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also skip points that have non-universal behavior when 
fitting. For our fits in this paper, all points are included. 

Another issue arises in regions where one has sparse 
data; there may be bins that do not have any obser- 
vations. For these zeros, the error bar should not be 
zero! We can use maximum likelihood methods to es- 
timate theoretical errors. Say we have N experiments, 
and bin the data with Li sizes in each bin i. With the 
median size in the bin Si, the expected probability for 
each size is pi = A(Si)/Si, where A(S) is the theoret- 
ical distribution of sizes, and so getting a zero in one 
of the bins will have the probability p = (1 — pi) NLi . 
Using maximum likelihood, this results in a cost of 
C = — \ogp = —NLilog(l — pi) ~ NLip given that p 
is small. The residual we add to the total cost is there- 
fore: 



Tzero = y/NLipi(0). 



(A4) 



This calculation is generalizable to sparse data that is 
non-zero- say one has n events in a bin for N mea- 
surements, then the probability is p = ^^L^-ny. — 
Pi) NLi ~ n Pi i an d the appropriate error bar we get is ap- 
proximately a = ty^t- One could then use the larger of 
the error bar given by this argument or the statistical 
error bars from the simulation. In practice, for our sim- 
ulations, we find it sufficient to use the statistical error 
bars but to compensate with minimum error bars given 
by Eq. \M 



In non-linear least squares fitting, we can also include 
priors in the cost if we have assumptions or information 
a priori about the parameters. For example, in our prob- 
lem, we put priors on the exponents inside the scaling 
functions (noo, ^oo? etc) to prevent them from going to 
large values and forcing their corresponding coefficients 
to zero. As a result our cost function now becomes: 



C(0) = -s* 



theory 



yf ata {9) 



S^ + log^) 



(A5) 

In Eq. |A5[ n represents all the arbitrary exponents that 



occur inside the scaling functions (n z 



m 



zzi $zz 



etc.). 



\og(Sf 1 ) is included to prevent Sn from going to zero. 

A clear understanding of these techniques for fitting is 
important for increasing the reliability of our results, and 
acknowledging its limitations. In a following subsection, 
section |A4[ we will further discuss how to estimate the 
reliability of results inferred by fitting data to a theory, 
generating systematic error bars for fitting results. 



3. Software for fitting 

To facilitate the exploration of this prob- 
lem, we have developed a software environ- 
ment, Sloppy Scaling, in Pytho n. This code is 
downloadable at http://www.lassp.cornell.edu/ 
sethna/Sloppy / SloppyScaling/SloppyScaling.htm. The 



main features of this code include various nonlinear-least 
squares fitting methods [38 , automatic plotting for 
visualization, and methods for generating systematic 
error bars on the theory. 



4. Systematic Error Bars 

We have quoted in our results systematic error bars 
instead of the more commonly used standard error bars in 
the parameters. Standard errors given by the covariance 
matrix are expected to be erroneous for our problem, 
since our problem is highly nonlinear in the parameters, 
and also sloppy- parameter combinations in the sloppiest 
directions can vary an infinite amount without affecting 
the fit. We use a method due to Frederiksen et. al [35] 
for Bayesian estimation of errors. This method involves 
assuming that given a theory (M) which is imperfect, 
a spread of parameters (each corresponding to different 
models) may fit the data(D) in an equally acceptable 
matter. We can define a probability of a certain model 
with 



P(6\D,M) =exp(-C(0)/T) 



(A6) 



where C(0) is the cost at a given set of parameters 0, and 
the effective temperature T sets a scale for the fluctua- 
tions away form the best fit. Since the cost at the best 
fit parameters Cbf is a measure of how well the theory is 
doing, we choose T = /N, where N is the number of 
parameters with "equipartition" allowing each degree of 
freedom a contribution of \T to the total cost. 

Ideally, after defining such a probability, one should 
sample parameter space to determine the systematic er- 
ror bars on parameters. However, in our functions sam- 
pling is non-trivial due to parameter evaporation [39] , 
the "entropy" for parameters drifting to infinity over- 
whelms their cost in degrading the resulting fits. There- 
fore, we make a quadratic estimate of the fluctuations in 
predictions, essentially using propagation of error to es- 
timate the systematic error. The covariance matrix gives 
an error a sta t that assumes the temperature of 1, corre- 
sponding to P(Q\D,M) = exp(— C{0)). Using propaga- 
tion of error, we calculate the systematic error according 
to our effective temperature T = 2Cbf/N: 



J sys 



Vfa c 



(A7) 



The shaded plots are generated by sampling according to 
the Hessian at the best fit, weighting our steps in each 
eigendirection by the inverse square root of the eigen- 
value, and scaling the steps with a low temperature (Tc). 
Then for our ensemble of parameters we calculate the 
fluctuations in the theory (residuals Sr ens ) correspond- 
ing to the ensemble. We scale up these fluctuations ac- 
cording to the temperature defined by the best fit, or 



Sri 



dr 
Z-50 T 



(A8) 
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We have also estimated systematic errors by removing 
large k curves (which we believe have larger corrections), 
and looking at the corresponding drift in exponents and 
parameters. The estimate of the systematic error that 
this procedure gives is often similar (sometimes smaller) 
to the one using our temperature-scaled propagation of 
error. 



5. Corrections to scaling 

Corrections to scaling play an important role in 
our scaling functions, their inclusion accounts for non- 
universal effects, and helps increase the reliability of uni- 
versal predictions. In each of the scaling functions in this 
paper we have included analytic corrections to scaling 
that capture the lattice effects of our automata simula- 
tions. They are of the form: 

exp(A7S + A 2 /S 2 ) (A9) 

This expansion for small S corrects for the lattice effects 
on small size avalanches which we believe to be present 
in the distributions. One can imagine that experiments 
may have other origins of non-universal effects, such as a 
nonlinear signal amplifier or distorting lens; one should 
always attempt to account for and include these [30] . 

One caveat is that, just as adding extra free parame- 
ters does not necessarily increase the quality of one's fits, 
adding corrections is not a guarantee for increasing the 
accuracy of one's scaling function. One should be careful 
in checking that the terms included in the corrections to 
scaling behave as expected in the region of interest, are 
subdominant when taking the appropriate limits, and do 
not confuse the main universal scaling function, either by 
canceling out terms or having the same effect. An exam- 
ple of this complication is seen in our studies of using the 
limit of the 11 distribution at W = for the local height 
distributions. 

In sections |IVD| and |VC[ we take the view point of 
using the parameterized from of An(s\Lk, W = 0) as the 
proper limit of A(h x \Lk). However, this data should be 
matched as well by Au(s\Lk, W = 1), since in effect this 
is what we are fitting for the local height distributions. 
We can view the ratio between the two functions as a 
multiplicative correction to scaling from lattice effects 
(without considering the corrections to scaling for lattice 
effects in equation |A9[ ): 

A n (s\L k ,W = l) 



f(h x , L k ) 



An(s\L k 
exp(-(Z7n( 



W = 0) 

h x N i/2 



(A10) 
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Here we define h s = h x /L k . Note however that for the 

term CnL^' 7111 ^ 2 h' 171 * , the value of the fit for n n 
equals 1.65, which is within the range of error for 1 + ( = 
1.63±0.02. Therefore this term is nearly Cn^ m2 , where 
m2 = 1.61 ± 0.6. This term then has the same effect as 
the correction to scaling term A^ 1 / s 2 at small W (since 
s ~ hW for 11 avalanches at small W). Since A^ 1 / s 2 is 
only significant in the range of small s, which only occur 
in the 11 distributions at small W, these two parameters 
serve the same purpose, and it is redundant to include 
both for the fits. We therefore remove the term A^/s 2 
for the 11 distributions for our fits. 

One can also check the corrections and see all pow- 
ers of L/c are negative, and the multiplicative correction 
approaches unity as we get closer to the critical point. 

— 1/2 

With our current fit, the term with L k dominates the 
corrections. Notice that, originally, to account for lat- 
tice effects, we have added corrections in integer powers 
of h x (ex.p(A\ 1 /h x + Al 1 /h 2 )), which are subdominant 
compared to L k X ^ . This implies that there are more 
dangerous corrections to scaling than originally inferred. 
By this method, one might check if the corrections orig- 
inally included are sufficient- in this case they are close, 
but could be improved upon by systematically adding 
similar terms. 

In a continuum case, the concept of local heights h x 
should describe a smooth shape tracing the depinning 
line, in our automata, since the smallest width is nat- 
urally 1 lattice spacing, the smallest possible window 
width is W = 1. In measuring the local heights, this 
discreteness limits the smoothness of the shape of h x , 
and gives rise to the corrections we see. Here we've seen 
that if we "know" from other measurements (in our case 
the 11 spanning avalanches) the right limit the universal 
scaling form should take, we may find the form of the 
corrections. 



6. Scaling collapses and their limitations 

We have argued at various points in the paper that 
scaling collapses are limited and may lead to questionable 
results. Here we will illustrate an example of this. Using 
the critical exponents given by a free fit of the windowed 
distributions (£ allowed to vary), we collapse the sizes, 
heights and widths. Comparing the figures included here 
and the ones in Section IV (Figures [5j [6j and [7[), one 
can see that the collapses are of similar quality. 

The fact that such distinct values of £ can yield sim- 
ilar quality collapses may imply (1) our "systematic er- 
ror" bars on £, estimated to be ±0.02, are in reality much 
larger, (2) collapses do not incorporate non-universal cor- 
rections to scaling, and these may have an important ef- 
fect, (3) collapses are not a reliable way of verifying the 
values of critical exponents. In particular, we expect cor- 
rections to scaling due to large k to be responsible for 
the drift in exponents. 
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FIG. 16: Size distribution collapse (color online) Here 
we collapse the size distributions with exponents r = 1.25, 
z/fc = 0.44, and £ = 0.68. Notice that the collapses are similar 
to the ones shown in Figure [5] Here only the combination of 
z/fc(l + C) affect the scaling collapse, the large shifts in Vk and 
f mostly cancel in the product, yielding similar collapses. 



FIG. 18: Width distribution collapse (color online) Here 
we collapse the size distributions with exponents r = 1.25, 
Vk — 0.44, and £ = 0.68. Notice that the collapses are similar 
to the ones shown in Figure [7] Here only the combination 
of z/fc affects the scaling collapse and (r — 2)(1 + Q affects 
the shape of the scaling collapse. Since Vk does not change 
significantly, the quality of the collapses are similar. 
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FIG. 17: Height distribution collapse (color online) Here 
we collapse the size distributions with exponents r = 1.25, 
z/fc = 0.44, and = 0.68. Here only the combination of h>kC 
affect the scaling collapse, the large shifts in and ( mostly 
cancel in the product, yielding similar collapses. Comparing 
this with the collapse shown in Figure [6] we see that the two 
collapses are comparable in quality, where in Figure [5] the 
large avalanche cutoff is collapsed nicely, and here the smaller 
avalanches are collapsed better. 



With these software tools and analytical methods, data 
at critical points may be analyzed while including mul- 
tiple scaling variables, allowing for the treatment of a 
broad range of experiments, and also allowing for a far 
more rigorous estimation of statistical and systematic er- 



rors. By using automatic fits to entire scaling functions, 
instead of traditional collapses, and by estimating sys- 
tematic error bars, we facilitate the interpretation of data 
with multiple scaling variables and analytic corrections 
to scaling. This advance will allow for better characteri- 
zation not only of noise in magnetic thin films and similar 
avalanche dynamics, but should be broadly applicable to 
all applications of critical phenomena and scaling theo- 
ries to experiments and simulations. 



Appendix B: The roughness exponent £ 

In the investigations presented in this manuscript, we 
have found that the estimates of the critical exponent £ 
have been problematic. In this appendix, we will discuss 
various means of measuring this exponent, the signifi- 
cance of the range of values we find from various measure- 
ment methods, possible origins of this range, and impli- 
cations for future research. We emphasize that any value 
of C in the range we observe (0.62 ± 0.02 to 0.72 ± 0.02) 
can describe all of our data essentially equivalently well. 

The shape of the front has been studied as an identify- 
ing feature for front propagation models, which is usually 
characterized by defining a roughness exponent which 
is measured through a height-height correlation function: 



C(r) = ((h(x + r)-h(x)) 2 )~r 2(: 



(Bl) 



The quenched KPZ model we use is conjectured to be- 
long to the Directed Percolation Depinning (DPD) uni- 
versality class [9], which is conjectured in turn to belong 
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FIG. 19: Height-height correlations for the qKPZ sim- 
ulations (color online) Shown here are the roughness ex- 
ponents for various simulation sizes L and k. We measure 
the height-height correlation function C(r) ~< (h(x + r) — 
h{r)) 2 >. A power law fit shows £ falls between 0.63 and 
0.68. The lower red-dashed line shows £ = 0.63 and upper 
black-dashed line £ = 0.68. The lines were shifted to show 
each individual power law. 



to the Directed Percolation(DP) universality class. For 
a pinned interface in DPD, the roughness exponent £ = 
0.63 ± 0.01 [lOj [11] matches that of the ratio of correla- 
tion length exponents in DP u±/u\\ = 0.6326±0.0002 [40] . 
(Here (,± ~ \p — p c \ u± and £y ~ |p — p c |^n , where p is the 
branching probability, and p c is the percolation thresh- 
old.). One may note that for a moving interface, the pic- 
ture is less well known; there has been numerical study 
that shows £ = 0.70 [EJ|41], but also arguments that the 
interface under this condition is not self-afflne [13], that 
the moving regions have ( = 1 and the pinned regions 
have £ = 0.63. One could imagine that our demagne- 
tizing force, like the velocity in DPD, could lead to a 
heterogeneous mixture of different scaling regions, con- 
verging to £ = 0.63 as k = 0. 



Figure 19 shows measurements of the height-height 
correlations in our model. Using finite-size scaling for 
a numerical fit, we see that as k is tuned away from zero, 
£ falls between 0.64 and 0.68, increasing with k. Measur- 
ing the local log-slope, one can see clearly a drift in the 
measured exponent in Figure [20} 

The range of our estimates (varying from £ = 0.63 ± 
0.02 to £ = 0.72 ±0.02) is large compared to our error es- 
timates; however, best fits with £ fixed within this range 
had costs within 1.5 times that of the best fit value, indi- 
cating both that our quadratic estimates for the system- 
atic errors are too small and that it may be challenging to 
definitively measure £ in either simulation or experiment. 

In our fits we find that C = 0.62 ± 0.02 for the 



size, width and height distributions joint fit, and ( = 
0.68 ± 0.02 for the 00, 10, 11 joint distributions. No- 
tice that although the direction of front propagation is 
in general along the y-axis in our problem, portions of 
the front will be at various angles to the y-axis. Since 
the local direction of the front propagation is not fixed, 
we can also choose a rotationally invariant definition of 
height and width: defining the height and width of an 
avalanche along the axes of the moment of inertia tensor. 
We define root-mean-square heights and widths as the 
square root of its eigenvalues. So if we fit rms heights 
and widths jointly, we get £ = 0.72 ± 0.02 much higher 
than 0.62 measured along the global axes. The difference 
in these two exponents seems to indicate that the local 
avalanche shape has a different geometry than the global 
avalanche front. Many front propagation models sponta- 
neously break rotational symmetry through the orienta- 
tion of the front. (Envision a circular front growing from 
a point, with differing front orientations.) Note, how- 
ever, that the qKPZ dynamics is anisotropic, breaking 
rotational invariance. 

We have also looked at fits with the windows scaling 
functions involving subsets of simulations with different 
k. Using smaller values of k generally lead to fits of £ 
closer to 0.63. This may point to imperfections in our 
theory function (are we missing some of the scaling be- 
havior dealing with fc?), or corrections to scaling (analytic 
or singular). One possibility is that there is a residual 
crossover effect having to do with another relevant vari- 
able. In the anisotropic form of the qKPZ model, the 
nonlinear term (Eq. [4]) A is a relevant variable, and is 
non-zero under renormalization [3j El [42] , and although 
we simulate the model at fixed A, we have observed there 
is a crossover effect following a direction having to do 
with both k and A. We note that without the nonlin- 
ear KPZ term the qKPZ model becomes the quenched 
Edwards- Wilkinson (qEW) model, and that for the qEW 
C > 1. Literature suggests there also may be a crossover 
effect due to a runaway fixed point [42] . We think these 
last two possibilities can be explored with more simula- 
tions done on different A and examining a crossover 
to the linear version of the qKPZ model (the quenched 
Edwards- Wilkinson) , to make a more complete picture 
of the phase space. 
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FIG. 20: Roughness exponents for the qKPZ simulations (color online) Shown here are the measurements of the local-log 
slope ln[C(r)/C(r/2)]/ln(2) of the height-height correlation function; this is a measure for the roughness exponents for various 



simulation sizes L and k. The results of the local-log slope ln[C(r)/C(r/2)]/ ln(2) is consistent with what is seen in Figure 19 



The lower red-dashed line shows f = 0.63, corresponding to directed-percolation depinning (DPD) and suggested by literature 
to be the correct value for our model. Whereas the upper black-dashed line shows ( = 0.68 which is the result of our fits 
of windowed avalanche distributions. In the blowup of the region of r/2 = 10 — 100, we can see there is a trend of larger £ 
corresponding to larger k simulations. 
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